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I. INTRODUCTION 



It has been known for thirty years that quantum field theory simplifies enormously if the 
number N of internal field components tends to infinity. In the case were the N components 
form a vector this leads to exact solutions in any dimension of space-time. For physical 
applications, ranging from solid state physics to gauge theories and quantum gravity, a 
different situation is much more pertinent: The case of iV 2 internal components that form 
a matrix. Here exact solutions have only been produced for very low dimensionalities. It 
is one of the outstanding problems of theoretical physics to extend large N technology to 
physically interesting dimensions. 

In the present article we will be concerned with matrix "spin systems", that is D- 
dimensional Euclidean lattice field theories whose internal degrees of freedom are hermitian, 
complex or unitary N x N matrices. The idea is to treat the problem by a three step 
procedure: 

(1) Eguchi-Kawai reduction: Replace the N = oo field theory by a one- matrix model coupled 
to appropriate constant external field matrices. 

(2) Character expansion: Express the partition function of the one-matrix model of (1) as 
a sum over polynomial representations - labelled by Young diagrams - of U(N). 

(3) Saddle point analysis: Find an effective Young diagram that dominates the partition 
sum of (2) in the large N limit. 

The insight that step (1) is possible is due to Eguchi and Kawai [|IJ. Intuitively it 
says that, if a saddle point configuration exists at iV = oo, it should be given by a single 
translationally invariant matrix (the so-called master field). In practice the reduction is 
rather subtle, and we will be using the twisted EK reduction which results in a one- 
matrix model in external constant fields encoding the original (discrete) space-time. 

Step (2) is novel in this context and is the main focus of the present work. The one- 
matrix model of (1) still has N 2 degrees of freedom, and it is well known that a saddle point 
for matrix models can only be found once the degrees of freedom are reduced as iV 2 — > N. 
The external fields encoding space-time prevent any naive reduction to the iV eigenvalues of 
the matrix, which is the route of choice for simpler models without external fields. But it is 
possible to replace the matrix integral by a sum over partitions corresponding to a sum over 
all polynomial representations of U(N). The crucial point is then that one ends up with 
a kind of one-dimensional spin model in Young diagram space with only N variables: the 
lenghts of the iV rows of the diagram. 

Step (3) might appear to be an exotic idea: we claim that the N = oo "master field" can 
be described by a "master partition" . However, it has already been recently demonstrated in 
a series of papers H, M that certain infinite sums over partitions are dominated by a saddle 
point configuration. This led to the solution of matrix models in external fields not treatable 
by any other method. The present models are more complicated, but not fundamentally 
different. 

The character expansions we find lead to a very interesting and apparently novel com- 
binatorial problem in Young pattern space (see section 4). More insight into this problem 
will be needed in order to proceed with the final step (3) of our program, the saddle point 
analysis. We introduce what we call "lattice polynomials" S/j,T/j which are polynomials in 
jf. They depend on the Young diagram h and the precise nature of the space-time lattice. 
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It might be objected that the present approach is futile unless one can demonstrate that 
the lattice polynomials H^,T^ can be explicitly computed or at least bootstraped at A = 
oo. But there is one important argument against this pessimistic assessment: The lattice 
polynomials H^,T^ only depend on the nature of the lattice but not on the local measure 
of the minimally coupled (matrix) spins of the modelQ. Therefore, solving interacting field 
theory in our language is of the same degree of complexity as solving the free field case. 

Finally we should mention that our program is very general since it applies in principle 
to any large N matrix spin system. It would be interesting to extend the method to matrix 
field theories with a gauge symmetry such as Yang-Mills theory. Indeed the EK reduction 
was initially designed for lattice gauge theory |TJ. Recently it was demonstrated by Monte 
Carlo methods that even the path integral of continuum gauge theory may be EK reduced 
to a convergent ordinary multiple matrix integral ||. A rigorous mathematical proof, as well 
as an investigation on whether the reduced model reinduces the field theory as N — > oo, are 
still lacking. At any rate, reducing a D- dimensional gauge theory, one so far ends up with 
a nonlinearly coupled D-matrix model, which is not yet tractable by the present machinery 
unless it is understood how to perform a further reduction DN 2 — > N 2 . 



II. REDUCED MATRIX SPIN SYSTEMS 



Consider a spin model on a periodic lattice. In order to be specific we will scetch the 
method for a two-dimensional lattice, but higher (or lower) dimensions can be treated as 
well. We will not dwell on details since they are well explained elsewhere. The variables are 
N x N hermitian matrices M(x) defined on the lattice sites x 



Zn= fl[VM(x) e 



S n = ATr £ 



-M(x) 2 + V(M(x)) - - J2 [M{x)M{x + A) + M{x)M{x - fi)} 
^ 2 M=12 



where fi denotes the unit vector in the /i-direction. It is equally natural to consider general 
complex matrices $(x) G GL(A, C), in which case 



GL 



x) e 



-<Sgl 



iSgl = YTr ]T + V($(x)& (x)) [$(x)&(x + jl) + $(x)$ f ( 



X 



(2) 



Except for the global symmetry of the matrix spins. In this paper we develop the theory in 
parallel for the case of U(iV) global symmetry (hermitian matrices) and U(Y) x U(iV) symmetry 
(complex matrices). The other classical groups could presumably be treated as well, but it is well 
known that they do not lead to different large N limits. 
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If V = in eqs. ([[]), (Q) the model is free. The integration measures in eqs.(P,(||) are the flat 
measures for hermitian and complex matrices: 



A? 



DM = ]J 



dMi. 



1 r dReMijdlmMij 



dRe$ ii dIm® ii 



ttN- 



(3) 



A third, very important type of spin model is the so-called chiral field, which looks like the 
free complex model eq.(fj) 



S v = -pNTi 



[U(x)tf(x + A) + U(x)U\x - p,)] . (4) 

but the matrices U(x) G U(iV) are unitary. In this case the measure T>U(x) is the Haar 
measure on the group. The model is therefore not free. 

The Eguchi-Kawai reduction , states that the above lattice models can be replaced 
at TV = oo by, respectively, the following one-matrix models coupled to constant external 
field matrices P and Q: 



Zu- 



DMexp iVTr 



-M 2 — V(M) + 0(MPMPi + MQMQ ] 



D$ exp iVTr [ - - V($$ t ) 



x 



exp j3NTi ($P$ t P t + $P+$tp + $g$tQt + $gt$tg 



Zu = /" Pf/exp flNTi (uPtfP* + UP^U ] P + UQU ] Q ] + UQ^U ] Q 



(5) 



(6) 
(7) 



Here P 



Pv and Q : 
/0 



Qjv are the famous N x N unitary "shift and clock" matrices 



P 



AT 



V 



\ 



o 1 



Q 



N 



\ 



N-2 



V 



N-l 
N 



(8) 



where u>n = exp and PnQn = ^nQnPn- To be more precise, the free energies as well as 
appropriate correlation functions (see 0) are identical to leading order in in the lattice 
field theory and the corresponding one-matrix model. The thermodynamic limit, that is 
a lattice of infinite extent, is approached when iV — > oo. We see that the structure of the 
lattice has been "hidden" in index space! It is natural to generalize the situation to a toroidal 
K x L lattice: 



P = P 



K 



In, 

K 



Q = Q L 



In, 

L 



(9) 



where N is chosen to be divisible by K and L. This allows to take the thermodynamic limit 
and the large N limit independently. If we put L = 1 (we can then equivalently omit Q 
altogether) the target space becomes a closed one-dimensional chain. 

We suspect that matrix models on arbitrary discrete target spaces can be EK reduced 
by appropriate external matrices, but this has not been worked out yet. 
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III. CHARACTER EXPANSIONS 



Now we turn to step (2) and rewrite the reduced hermitian, complex and unitary matrix 
integrals eqs.@,(^),(0) as sums over representations of U(iV). To this end introduce the 
following source integrals: 



MA 



VM exp iVTr [ - -M 2 - V(M) + JM 



(10) 



Z GL [JJ] = / £>$ exp NTi [ - $$ f - F($$ t ) + J$ + $ f J 



(11) 



Z V [JJ] = VU exp iVTr \jU + iflj 



(12) 



The two different ways of introducing a source are due to the U(iV) symmetry of hermitian 
matrices on the one hand and the U(iV) x U(iV) symmetry of complex (and complex unitary) 
matrices on the other. The reduced models are easily obtained from the source integrals by 
applying an operator: 



Z u = exp £-Tr (dPdP^ + 8QdQ ] ) ■ Z U [J] 



N 



j=o' 



(13) 



Z GL ,v = exp -^Tr (dPBP^ + dP^BP + dQOQ^ + OQ^Bq) ■ Z GL)U [JJ] 



j=j=o 



(14) 



Here d,B denote N x N matrix differential operators whose matrix elements are dji = -^j- 

and = -t—-. It is clear that the source integrals are class functions of, respectively, J and 

J J. Therefore they can be expressed as character expansions, with known (see ||, ||, ||) 
expansion coefficients. If V — 0, they read for the hermitian and complex source integrals, 
respectively, 



Z H [J]=exp lNTrJ 2 = Y,Xh(A 2 ) X h(J), 
1 h 



(15) 



Z GL [JJ] = exp iVTr J J = £ Xh{M)Xh(J?), 



(16) 



while for the unitary source integral one has || 



h XhKV 



(17) 



Here the sum runs over all partitions h labeled by the shifted weigths hi = N — i + rrii, where 
Tn% > 0, i — 1, . . . , N, is the number of boxes in the z'th row of the Young pattern associated 
to h. Xh{J) is the Schur function, dependent on J, on the diagram h. It is identical to the 
Weyl character of the matrix J corresponding to the representation labeled by h. A± and 
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A 2 are defined through Tr A\ = N(5 k ,o + and TrA^ = N(5 k ,o + ^,2), and XhW is the 
dimension of the representation. For more details on the notation, and for explicit formulas 
for the characters Xh{A\), Xh(^-2) and XhiX) see 0) @- F° r a non-zero potential V, the 
hermitian and complex character expansions become a bit more complicated, but are still 
available: 



Zn[J}=J2 & hXh(J), 



(18) 



where 6^ is given by 



e. 



and flh by 



XhW 



, XhW 



z GL [jj] = J2^hXh(JJ), 

h 



(19) 



DM exp iVTr 



-M 2 -V(M) Xh(M) 



(20) 



£>$ exp TVTr f - $$ f - y($$ t )l 



(21) 



The integrals appearing in eqs.(|20| ) , (|2l|) are ordinary one- matrix integrals which may be 
computed rather explicitly as N x N determinants. Their analysis in the iV — > 00 limit 
proceeds by employing standard techniques, supplemented by the methods of ||. 

Now we apply the operators in eqs.(|l3|),(|T4l) in order to generate the space-time lattice; 
this results in character expansions for the reduced matrix field theories. In the hermitian 
case one has (here \h\ = J2i m i =number of boxes in the Young diagram) 



Z n = T,Xh(A 2 )E h (3~ for V = 0, 

h 



(22) 



with 



Z H = '£Q h Z h p* for V^O, 

h 



E h = exp iTr (dPdP^ + 8QdQ^) ■ Xh (J) 



j=o 



The free complex, interacting complex, and the unitary case become 

Z GL = Y,Xh(A 1 )T h ^ for V = 0, 

h 



(23) 



(24) 



(25) 



Zgl — J2^hT h p\ h \ for V^O, 



(26) 
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^ — win — fc ^ ' 



with 



T fc = exp ^Tr (dPdP* + dP^dP + dQdQ^ + dQ ] dQ) ■ Xh{JJ) 



j=j=o 



(27) 



(28) 



The character expansions eqs. (p2]) , (p3|) , fl2~5|) , fl26|) , fl27|) are at the heart of our proposal. It 
is seen that they neatly separate the nature of the local spin weight [xhi^-i) i^hiXh{^i) P-hi 
(X/i(^4i)) 2 (Xft(l)) _1 ) an d the nature of the embedding space (H/,,T/,). As a striking example, 
note that from the point of view of our character expansion method the difference between 
the free Gaussian model on a toroidal lattice eq.fl2"5|) and the non-trivial chiral model eq.(f2~7|) 
is a simple, explicitly known factor 



Xh(A 1 



N^ N 

i=i 



hA 



The character expansions involve sums over TV variables only and we can write down a saddle 
point equation for the effective density of the master partition. In order to complete the 
program, we need a second bootstrap equation for the novel quantities E h and T^, which 
contain the connectivity information of the lattice. 



IV. LATTICE POLYNOMIALS 



Inspection of the quantities E h and T h in eqs. (|2"4]) , (f2"8|) shows that they are polynomials 
in the variable of degree not higher than, respectively, — 1 and \h\ — 2. They are zero 
if the number \h\ of boxes in the Young pattern is odd. Conjugating the diagram gives the 
same polynomial except for the replacement j? — > — The first few can be computed by 
brute force calculation directly from the definitions eqs.fl24]),(|28|), see Table. 



Table : The first few D = 2 lattice polynomials 



h 




r h 


2 


2 


2 


l 2 


2 


2 


4 

31 

2 2 


3 + 12^ 

5 + 4^ 
6 


3 + 24^ + 54^ 
5 + 8^ + 18^ 
6 


21 2 
l 4 


5-4^ 
3-12± 


5-8^ + 18^ 
3 - 24^ + 54^ 



Here we used Tr (P k Q l ) = NSkflSifi, which is true as long as \k\ < N, \l\ < N. We also 
replaced to^ — > 1, uj* n — > 1 (remember uj n = exp ^p): in other words, we assumed P and 
Q to commute at large N. Both assumptions are innocent at least in the strong coupling 
(small pi) phase. If the model possesses a weak coupling phase (like e.g. the chiral field 
eq.(|7|)), these assumptions may have to be reconsidered, if we want the character expansion 



7 



to describe this second phase as well. This is because in the present approach we expect 
large iV phase transitions to correspond to the situation where the number of rows of the 
master partition is of O(N) ("touching transition"). Note that we cannot drop the other 
terms of O(j^) in S/j,T^ since the character expansions are for the partition function and 
not for the free energy. 

The direct calculation of the lattice polynomials quickly gets very tedious. The com- 
binatorics involved seems to be of a novel type. While we have not yet found an efficient 
calculational scheme or recursive method, let us give some interesting representations for H/, 
and T h that may prove useful later. Introduce the following Gaussian measure on the space 
of M x iV (M < N) complex matrices A: 

m = g g ( dReAjjdlmAjj ^ ^ ^ ^ _ ^ (2Q) 



i=l 



j=l \ 



This measure is invariant under U(M) x U(N). It is then fairly easy to prove (c/ |Q) the 
following representation for the character of the source: 

Xh (J) = J 'vU Xh (tf) J [DA] ex P NTrUAJA\ (30) 

where U G U(M) is unitary and T>U is the Haar measure on U(M). This formula is valid 
for diagrams h with at most M rows. Therefore becomes, cf eq. (|2H]) 

E h = J VU X h{U^) J [DA] expNTr ^UAPA^UAP^ + A ] UAQA ] UAQ^. (31) 

After a Hubbard-Stratanovich transformation decoupling the quartic terms by Gaussian 
MxM complex matrices S and T (with measure as in eq.(^) with A" — » M), and integration 
over A, we obtain the representation 



~ h = j VUxhiU^) j[VS][VT]x 

00 1 _. ....... ...... 

Jfc 



x exp 



Tr M ^ v £-(S[/®P + S' t £/®P t + Tfy®Q + T t [/®Q t ) • ( 32 ) 



k=l 



The combinatorial interpretation of the exponential in eq.(^) is the following: we have a 
generating function for a non- commutative random walk on a two-dimensional lattice with 
variable U. The representation is useful for getting some exact results on the H^, but we 
have not yet been able to compute the integral eq. (|32"D exactly except for M = 1 (characters 
with just one row). E.g. we can find a generating function (with Z{ being the eigenvalues of 
U) for the large A" limit of 



M i 



ZiZj) 2 



E^ = °°X^) (33) 



giving the constant terms of the lattice polynomials. This is however not sufficient for the 
large A" limit of the field theory, as already mentioned. A curious feature of eq.(^) is that 
we can take A" — > oo while keeping M in the range 1 < M < JV. That is, it should be 
possible to find a saddle point for the situation where the row lenghts are large compared 



S 



to the number of rows, corresponding to the extreme strong coupling limit. Furthermore, it 
should be investigated whether the M x M matrices can be taken to commute as A" — > oo. 

Similar, if slightly more complicated representations are possible for T^; here the starting 
point is the expression 



Xh(JJ) = J VU X h{U ] ) J [PAi][PA 2 ] expiVTr (XJ^A X JA\ + A 2 JA\U^ 
which means the lattice polynomials become 
T h = f VU X h(tf) j[VA x )[VA 2 ]x 

X 



(34) 



exp NTi (a\U^ A X P A\U^ A 2 P^ + A\U^ A X P^ A\U^ A 2 P 
exp NTi (a\U^ A X QA\U^ A 2 Q^ + A\U? A X Q^ A\U^ A 2 Q 



(35) 



and the non- commutative random walk representation is 
T h = J VU Xh (U j ) J [VS][VS][VT][VT]x 

OO -y 



x exp 
x exp 



x 



Tr m ®n E T.( SU ^ ®P + SU l i(g)P^ + TU l 2(g)Q + TU* <g> Q 
k=i fc 

OO 1 , ■ 

Tr Mm £ - (3*11* g, P + S-t[/| g, F t + T+C/i ®Q + ® gt) 

fc=i 



(36) 



from which we find that 00 



^7V=oo 



cf eq.(|33|), but corrections are different (see 
Table above). Again, for arbitrary one-row representations (M = 1) it is possible to obtain 
Th rather explicitly. 

Another potentially useful representation^ of the lattice polynomials is given by the 
following dual equations: eq.(El) becomes 



Xh (d) ■ exp llr (jPJPt + jqjq 



J=0' 



and eq.(EH) is dual to 



T, 



Xh(dd) ■ exp — Tr (jPJP^ + JPUP + JQJQ^ + JQUQ 



j=j=o 



(37) 



(38) 



We could go on and discuss correlation functions which are naturally included into the 
present formalism. In particular, it is straightforward to give expressions for their character 
expansions in terms of modified lattice polynomials, and it remains true that the combi- 
natorics is independent on whether the reduced field theory is free or interacting. This is 
however beyond the scope of the present article. 

While it is unclear whether the D > 2 lattice polynomials can be computed exactly for 
a general partition, it should be stressed once more that this is unnecessary; all we need is 
an indirect method in order to extract the large A" behavior. 



2 We thank D. Verma for pointing this out to us. 
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V. CONCLUSIONS 



This solution to the problem of the large N limit of (non-gauge) matrix field theories is 
not yet complete since the structure of the lattice polynomials we introduced still needs to be 
further analyzed in order to be able to write the full set of saddle point equations. However 
we feel that we are definitely closing in on the large iV problem, and that we have brought 
it into the simplest form to date. The proposed approach is concrete, systematic and rather 
general: we demonstrated that the reduction from N 2 to N variables is possible once one 
changes variables from matrices to partitions. In this language the master field becomes a 
master partition. Presumably one should first (re) derive in the current framework the exact 
solutions for some lower dimensional target spaces before dealing with the two (and higher) 
dimensional field theories. 
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